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Abstract 

Phenomenological AdS/QCD models, like hard wall and soft wall, provide hadronic mass spectra 
in reasonable consistency with experimental and (or) lattice results. These simple models are 
inspired in the AdS/CFT correspondence and assume that gauge/ gravity duality holds in a scenario 
where conformal invariance is broken through the introduction of an energy scale. 

Another important property of hadrons: the decay constant, can also be obtained from these 
models. However, a consistent formulation of an AdS/QCD model that reproduces the observed 
behavior of decay constants of vector meson excited states is still lacking. In particular: for radially 
excited states of heavy vector mesons, the experimental data lead to decay constants that decrease 
with the radial excitation level. 

We show here that a modified framework of soft wall AdS/QCD involving an additional dimen- 
sionfull parameter, associated with an ultraviolet energy scale, provides decay constants decreasing 
with radial excitation level. In this version of the soft wall model the two point function of gauge 
theory operators is calculated at a finite position of the anti-de Sitter space radial coordinate. 
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I. INTRODUCTION 


The discovery of the AdS/CFT correspondence 1-3] - an exact duality between string the¬ 
ory in certain higher dimensional geometries and supersymmetric SU(N) conformal gauge 
theories with large number of colors - motivated the development of the AdS/QCD phe¬ 
nomenological models. One of the ideas behind the AdS/QCD approach is to assume that 
there is an approximate duality between a field theory living in an anti-de Sitter background 
deformed by the introduction of a dimension-full parameter and a gauge theory where the 
parameter plays the role of an energy scale. 

The idea of associating a cut off in anti-de Sitter space with an infrared energy scale 
in the gauge theory appeared in 4], where the AdS/CFT duality between a scalar field in 
anti-de Sitter space and a scalar operator (glueball) in a gauge theory was supplemented 
with the introduction of a cut off in the radial AdS coordinate. The amplitudes for fixed 
angle scattering calculated in this way show a “hard” behavior - they decrease with a power 
of the energy - as it is observed experimentally for hadrons and is obtained in QCD m 
One remarkable outcome of ref. |4j, was to show that, contrarily to what happens in flat 
space, where the amplitudes decrease exponentially with the energy, string theory in anti-de 
Sitter space is consistent with the physics of strong interactions, concerning the large energy 
behavior of fixed angle scattering amplitudes. 

This idea of placing a geometrical cut off in AdS space was then used in [7, s| as a model 
for calculating the masses of glueballs. The same approach was then extended to other 
hadrons and named as hard wall AdS/QCD model (for a recent review and a wide list of 
related references, see q| ). 

The hard wall model fits the available spectrum of radial excitations predicted by lattice 
for scalar glueballs [ 7 , 8] and experimentally observed for light mesons and baryons 10 ], 
However, the asymptotic behavior of the spectrum for the hard wall is that the masses grow 
linearly with the excitation number n for large n. In contrast, the available experimental 
data indicate that there is an approximate linear relation between mass squared and radial 
excitation number. 

Another AdS/QCD model, the soft wall, where the square of the mass grow linearly 
with the radial excitation number was introduced in [ll]. In this case, the background 
involves AdS space and a scalar field that acts effectively as a smooth infrared cut off. Many 
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interesting improvements of this model appeared then in the literature. In particular, for 
fermions the soft wall, as formulated in ref. [ll], does not lead to a discrete spectrum 
because the dilaton introduced in the action factorizes out in the equations of motion. This 


problem can be overcome, as shown in 
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14], considering that the mass of the fermionic 


held depends on the radial coordinate of AdS space. This change in the model is interpreted 
as effectively capturing the dependence of the dimension of the gauge theory operator on 
the energy scale, as it happens in a non conformal theory. Another important improvement 
was the discovery that using light-front coordinates, the hth coordinate of AdS space can 
be interpreted as representing the impact parameter jl5, 16]. Along the same line, the 
equivalence between QCD light-front and the equations of both hard and soft wall models 


was show in ref. 


17]. It is interesting to mention also that the soft wall model was extended 




to heavy quark hadrons in m 

AdS/QCD models provide also a tool for calculating another important property of 
hadrons: the decay constant. The decay of mesons, that in contrast to baryons are not 
subjected to a number conservation condition, is represented as a transition from the ini¬ 
tial state to the hadronic vacuum. The amplitude for the decay process is thus calculated 
from the holographic dual of the transition matrix between the vacuum and the one particle 
state. For the soft wall model, expressing the two point correlator as a sum over transition 


matrix elements, one finds the decay constants [Hi] (for a detailed discussion see 19)). The 
decay constant f n of a meson at radial excitation level n is related to the transition ma¬ 
trix between the one particle state and the vacuum. We use in this article the definition: 
( 0 | J M ( 0 ) | n) = e^f n m n . 

The results of the soft wall are degenerate: all the decay constants of the radial excitations 
of a vector meson are equal. For the hard wall model the decay constants of radial excitations 
are not degenerate but they increase with the excitation level. In contrast, the experimental 
results available for heavy vector mesons show that higher excited radial states have smaller 
decay constants. For the excited states of the light vector meson p the available experimental 
data does not provide an estimate for the decay constants. 

For the hard wall model the description of the decay process of vector mesons, as done 
in ref. 20] leads to f n oc 1 /J 1 ( 7 n ) where J\ is the Bessel function of order one and 7 n are 


the zeroes of the Bessel function of order zero: Jo (bn) = 0. These decay constants increase 
with the excitation level since the absolute value of J\ evaluated at the zeroes y n of Jq is 
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a monotonically decreasing function of n. So, the hard wall model does not reproduce the 
observed behavior of the decay constants. 

In the original soft wall model there is just one dimensionfull parameter, that plays the 
role of a mass scale. This infrared parameter is introduced in the dilaton background of 
the model and determines the mass spectrum. The main motivation of this letter is to 
propose that the description of decay processes in AdS/QCD should involve an additional 
dimensionfull parameter, or energy scale. In the usual AdS/QCD approach one describes the 
hadronic decays by means of held correlators calculated at the AdS boundary. Our proposal 
is to find the decay constants from the two point correlator calculated at some finite value 
z = z 0 of the radial coordinate of AdS space. This corresponds to introducing a new energy 
scale: 1 /z 0 in the model. 

A consistent way of calculating held correlators of the gauge theory operators at a finite 
position z = Zq of the radial AdS coordinate was presented in refs. [21,|22|. One interesting 
thing about this procedure is that the introduction of a second energy scale improves the 
description of the mass spectra of vector mesons made of heavy quarks. 

We specialize our discussion of the soft wall to the case of interest: vector mesons. The 
soft wall action for a vector held V m = (V/, 14) (/i = 0,1, 2, 3), assumed to be dual to the 
gauge theory current reads: 

1 


I = / d A xdz e ^ 


4 gi 


F F r ‘ 

2 mn 1 - 


(i) 


where F mn = d m V n — d n V m and <3> = k 2 z 2 is the soft wall dilaton background, that plays the 
role of a smooth infrared cut oh and A; is a constant representing the mass scale. 

The background geometry is anti-de Sitter AdS$ space, with the metric: 

ds 2 = e 2A ^ z \-dt 2 + dx-dx-\-dz 2 ), (2) 


where A(z ) = — log[z/R ) and the Poincare chart of AdS corresponds to (t, x) G 1Z 1,3 , z G 
(0, oo) and z is called radial coordinate. 

We use the gauge V z = 0. The boundary value of the remaining components of the 
vector held: V°(x) = lim z ^, 0 V^x, z), )jl = 0,1,2,3 are assumed to be, as in the AdS/CFT 
correspondence, the sources of the correlation functions of the boundary current operator 
J>*(x) (= q^q{;x)). Namely: 

(0| Jp{x)J v {y) |0) = fiyoij^ §yOv(y^ eX P ( _ lonsheU) > (3) 
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where the on shell action is given by the boundary term: 


L-on shell 



->—k 2 z 2 


-V tl d z V 11 


( 4 ) 


We have introduced gl = g\fR that is the relevant dimensionless coupling of the bulk vector 
held. For comparison with the proposal that we will present in the next section, let us 
emphasize that this action is calculated at the AdS boundary z — 0. One can then write 
the on shell action in momentum space and decompose the held as 


Vpfaz) = v(p,z)V%(p ), 


( 5 ) 


where v(p, z ) is the bulk to boundary propagator that satishes the equation of motion: 

_ k 2 z 2 2 

d z (— — d z v{p,z))+^e~ k2z2 v(p,z) = 0. (6) 

z z 

In order that the factor V®(p), defined in the decomposition of eq. (J5]), works as the source 
of the correlators of gauge theory currents, one must impose the boundary condition: 
v(p, z = 0) = 1. On the other hand, the two point function in momentum space is related 
to the current-current correlator by 

(p 2 Vuv ~ PuPv) n (p 2 ) = Jd* X e-^(0\M X )M0)\0). (7) 


Thus the two point function is holographicaly expressed in terms of the bulk to boundary 
propagator as: 


n (p 2 ) = 


95 ( -V 2 ) 


' v(p, z)d z v(p, z) 


( 8 ) 


The two point function has the following spectral decomposition, in terms of masses m 
and decay constants f n of the states of the gauge theory: 

oo r‘2 

J n 


n(p 2 ) = E 


5 (-p 2 ) - ml + it 


( 9 ) 


In the soft wall model the states are described by the normalizable solutions (f> n (z) of the 
equation of motion ([6]). One can represent the bulk to boundary propagator v(p, z) in terms 
of the normalizable solutions by introducing the Green's function 


G(p, z, z') 


hi ( -p 2 ) - m i ’ 


( 10 ) 
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and writing: 


v(p,z ) = 


e k2z ' 2 d z 'G(p : z', z ) 


( 11 ) 


The normalized solutions have the form: <j> n (z) = ^2(n — 1)! /n\(kz) 2 L) l _ 1 (A: 2 z 2 ) , where 
Ln(x) are the generalized Laguerre polynomials of order one. They satisfy the boundary 
condition lim z ^ 0 <f>n(z) = 0 and correspond to solutions of eq. (J6]) with momenta: — p 2 n = 
= 4 k 2 n , n = 1 , 2 ,.... 

The bulk to boundary propagator defined in eq. (fill) satisfies the equation of motion © 
and also the boundary condition v(p, 0) = 1. It is important to discuss the subtle way in 
which this boundary condition is satisfied. If one simply takes the definition for the bulk to 
boundary propagator from the Green’s function one finds an infinite sum of terms with the 
propagator poles l/(— p 2 — m 2 ) that would in principle lead to a dependence of v(p, 0) on the 
momentum p. However, the boundary condition liuq^o <j>n{z) — 0 leads to a sum of singular 


terms with vanis’ 
presented in ref. 


ring coefficients, that needs to be regularized. Such a regularization is 


19], where an integral representation for the bulk to boundary propagator 


is found. This integral representation is well defined in the interval 0 < z. It is explicitely 
equal to one in z — 0 and coincides with the expression (fill) when z > 0. 

The two point function obtained holographically in eq. dHJ) has the same pole structure 
of eq. (jU]) , with the residues corresponding to the decay constants. Thus: 


fn = 


d z <f) n (z) 


m n g 5 0 


z—>0 


ky/2 

95 


( 12 ) 


This result shows that in the soft wall model all the radial excitations of a vector meson 
have the same decay constant. This degeneracy contrasts with the decrease with excitation 
level n in the decay constants obtained from experimental data for heavy vector mesons. 
This problem was discussed in reference 23], where a finite temperature AdS/QCD model 


for quarkonium in a plasma was discussed. In order to describe the finite temperature 
behavior, the authors proposed a procedure to fit the zero temperature parameters. This 
was done by adding extra terms in the equation of motion for the normalizable modes, in 
such a way that one arrives at a model with four parameters that are fixed by the four 
quantities available: the masses and decay constants of J/ip and ?//. 
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II. ALTERNATIVE DESCRIPTION OF THE DECAY PROCESS 


The soft wall model, as reviewed in the previous section, contains one dimensionfull 
parameter: k, introduced in the dilaton background. This parameter plays the role of an 
infrared (or mass) scale of the model and determines the spectrum of the vector mesons: 
m n = Ak 2 n as well as the decay constants. The linear relation between the square of the 
mass and the excitation level is the expected asymptotic behavior. However, the observed 
behavior of decay constants of vector mesons is that they decrease with the excitation level 
while the soft wall model predicts degenerate values. 

The relation between the vacuum expectation value of the product of two gauge currents 
and the decay constants 

oo r2 

/ d 4 x e" 4p ' x (0| J^x) J u (0) |0) = (p 2 V^-p^Pu) J2 7—V— n 2 , • , (13) 

J V y ,7^i (-P 2 ) - rn 2 n + xe 


comes from introducing on the lefthand side of this equation a basis of the states of the 
theory. The decay constants are related to the transition matrices between the one particle 
states and the vacuum as: (0| J M (0) | n) = e^fnirin. Our point of view is that, in order to 
incorporate in the model the effect of the interactions that govern the transition processes 
from the one particle state to the vacuum, corresponding to the hadron decay, one needs 
to introduce one additional dimensionfull parameter in the soft wall. This can be done by 
calculating holographicaly the operator product of currents of eq. (TT3T) at a finite location 
z = z 0 of the radial coordinate. The new parameter, namely 1/zq, corresponds to an energy 
scale where the transition matrices are calculated. 

The idea of introducing an ultraviolet cutoff in the soft wall model was previously consid¬ 
ered in ref. 24], In this reference normalized solutions (j) n satisfying boundary condition at 
a position corresponding, in our notation, to a minimum value of the z coordinate z — 1/A 
were studied numerically. The condition imposed on the solutions was lim^i/A <j> n {z) ~ z 2 
; lim z _n /A <f>' n (z) rs j 2 z . The spectrum of masses obtained for the p meson from this ap¬ 
proach showed a better fit to the experimental data, in comparison with the original soft 
wall model. 

In this reference 24] the decay constants were also discussed. It was assumed that with 
the introduction of the UV cutoff they should be given (using our choice of radial coordinate 
and our convention for decay constants) by: 
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fn = 


d z 4>n(z) 


(14) 


z—tl/A 


m n g 5 z 

Using this relation, a decrease in the decay constants with the excitation level was found. 
However there is a problem with expression (1141) . The boundary value of the bulk held plays 
the role of the source of the gauge current correlators. When one takes the boundary to 
be at some finite value of the coordinate z — 1/A and decompose the held, as in eq. (J5j), 
Vfj,(p,z) = v(p, z)V®(p) , the source is V®(p) = linu^/A V/(p, z) and the bulk to boundary 
propagator has to satisfy the condition: lin^^/A v(p, z) = 1. In order to represent v(p,z ) 
in terms of the normalized solutions </> n in the same way as in equations (fTOj) and (HT| it is 
necessary that the condition lim,_ i . 1 /A 4>n(z) = 0, holds for the normalized solutions. This 
way, one finds for the bulk to boundary propagator at z = 1/A a sum of singular terms with 
vanishing coefficients that can be regularized to unity, as discussed in section II. 

For a non vanishing boundary condition, like the one used in ref. 24], the bulk to 
boundary propagator will be a sum of propagator pole terms with finite coefficients, keeping 
a momentum dependence. So, the condition: v(p, 1/A) = 1 is not satished and V®(p) is not 
the source of the correlators. Thus, in this case, one can not assume the expression (fT4l) to 
represent the decay constants. 

What comes out as a consistent approach to describe the decay constants by calculating 
the held correlators at a finite position z = z 0 of the radial coordinate is to follow a procedure 

I i n 

presented in refs. [21,[22j. In these references a cut off in the radial AdS coordinate was 
introduced so as to play the role of an ultraviolet cut off in the gauge theory. We will relate 
the same sort of cut off with the decay process. The bulk to boundary propagator is written 
as a solution of the equation of motion (J6]) but now excluding the region 0 < z < Zq. The 
solution is divided by a constant (in the 0 coordinate) so as to satisfy by construction the 
boundary condition v(p,z 0 ) = 1. Namely: 

U(p 2 /Ak 2 : 0, k 2 z 2 ^ 


^ U(p 2 /Ak 2 , 0, k 2 Zo) ’ 

where U(a,b,c ) is the Tricomi function. Using the new on shell action given by: 

1 r p -k 2 z 2 

T ___ / 

1 on shell l w 




-v,,d z v> 1 


>z 0 


we calculate the two point function: 

n (p 2 ) = 




v(v, z)d z v(p, z) 


9s(-P 2 ) 


(15) 


(16) 


(17) 


0 














Then, using the new bulk to boundary propagator of eq. (fl5l) we get: 


1 e k2z °U(l + p 2 /4/c 2 ,1, k 2 zl) 

2l 


(18) 


where a recursion relation for the Tricomi functions has been used. 

The spectrum of masses and decay constants comes from the analysis of the poles in the 
momentum variable p 2 . The Tricomi function in the numerator is non singular in p 2 , so the 
two point function of eq. (1T5|) is only singular at the zeros of the denominator. Furthermore, 
the poles p 2 of the function fl(p 2 ) are simple. So, in a neighborhood of p 2 = p 2 n one can 
approximate the two point function by: 


lim fl(p 2 ) 

p 2 ->p2 


f 2 

J n 


(-p 2 ) + pI 


(19) 


We associate the coefficients of the approximate expansion near the pole with the decay 
constant f n in analogy with the exact expansion shown in eq. ([9]) . This way we get the masses 
from the localization of the poles of the two point function and the decay constants from 
the corresponding coefficient. More precisely, let Xn be the roots of the Tricomi function: 


U(Xn,0,k 2 z$) = 0 , 


( 20 ) 


then the holographic vector meson masses are: m 2 = Ak 2 \n ■ The decay constants are 
calculated numerically from the fit to the approximate form presented in equation (fT9lh 
That means: 


fn = ] im 2 ((-p 2 ) +p 2 n ) n (p 2 


>p. 


( 21 ) 


The coupling g 5 = g 5 /\J i? of the vector field in the AdS bulk can be obtained by comparison 
with QCD, as explained in references [ll,[l9], wich gives: g§ = 2n. 


III. CONFRONTING THE MODEL WITH THE AVAILABLE DATA 


The experimental values for the masses m v and the electron positron decay widths 
T y_^ e + e — for vector mesons are taken from ref. 25]. The decay constant of a vector meson 


state is related to the corresponding mass and width by 

3tn v Ty >. e + e - 


26|: 


/ = 
J V 


47ra 2 c, 


(22) 
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Charmonium data 


Masses (MeV) 

T y _> e+e - (keV) 

Decay constants (MeV) 

15 

3096.916 ±0.011 

5.547 ±0.14 

416.2 ±5.3 

25 

3686.109 ±0.012 

2.359 ±0.04 

296.1 ± 2.5 

35 

4039 ± 1 

0.86 ±0.07 

187.1 ± 7.6 

45 

4421 ± 4 

0.58 ±0.07 

160.8 ±9.7 


TABLE I. Experimental masses and electron-positron widths from 


251 ] and the corresponding decay 


constants for the Charmonium S-wave resonances. 


Bottomonium data 


Masses (MeV) 

IV->e+e- (keV) 

Decay constants (MeV) 

15 

9460.3 ± 0.26 

1.340 ±0.018 

715.0 ±2.4 

25 

10023.26 ±0.32 

0.612 ±0.011 

497.4 ± 2.2 

35 

10355.2 ±0.5 

0.443 ± 0.008 

430.1 ± 1.9 

45 

10579.4 ± 1.2 

0.272 ± 0.029 

340.7 ±9.1 


TABLE II. Experimental masses and electron-positron widths from 
decay constants for the Bottomonium S-wave resonances. 


251 ] and the corresponding 


where a = 1/137 and c v for the two families that we consider are cj/^ = 4/9 and ex = 1/9. 
The uncertainties in the decay constants are: 

3 

5f v = 8 ?m 2 c f (m v sr v ^ e+e - +r v ^ e+e -5m v ) . (23) 

In table I we show the experimental values of the masses and widths and the associated 
decay constants with the corresponding uncertainties for the charmonium vector meson 
J/fp, made of a charmed quark anti-quark pair and for the first three radially excited S- 
wave resonances. This table shows that the decay constants decrease monotonically with 
the radial excitation level. 

Then, in table II we show the experimental data and the associated decay constants with 
the corresponding uncertainties for the botomonium vector meson T, made of a bottom 
quark anti-quark pair and for the first three radially excited S-wave resonances. As it 
happens for the charmonium, the decay constants decrease monotonically with the radial 
excitation level. 
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Now we present the results for the masses, obtained from equations (1201) and the decay 
constants obtained from equation (l2Tj) for charmonium and bottomonium S-wave states. A 
nice fit to the experimental data is obtained using the parameters 

k c = 1 . 2 GeV- k b = 3AGeV\ 1 /z 0 = 12.5 GeV, (24) 


where k c and k b are the values of the constants k used for charmonium and bottomonium, 
respectively. The values are different, reflecting the difference in the masses of charm and 
bottom quarks. The energy scale 1/zo is taken to be flavor independent and represents a 
characteristic energy of the decay process. In the non hadronic decay a very heavy vector 


meson annihilates into light leptons. So, 
that we consider. We show in tables III 


Results for Charmonium states 


Masses 

Decay constants 

15 

2410 

258.8 

25 

3409 

251.7 

35 

4174 

245.9 

45 

4819 

241.0 


TABLE III. Masses and decay constants in 
MeV obtained for the Charmonium states. 


the ultraviolet scale is of order of the heavier masses 
IV the results obtained for the charmonium and 


Results for Bottomonium states 


Masses 

Decay constants 

15 

7011 

627 

25 

9883 

574 

35 

12077 

538 

45 

13923 

512 


TABLE IV. Masses and decay constants 
obtained for the Bottomonium S-wave res¬ 
onances. 


bottomonium S-wave states, respectively. One can see that the model produces a decrease 
of the decay constants with the radial excitation level, as it happens with the available data 
for these particles. 

One interesting way to characterize the model in terms of predictability is to define the 
rms error for estimating N quantities using a model with N p parameters as: 




\ (N-N r )^\0. 


i {«>, 


(25) 


where Oj is the average experimental value and 50i is the deviation of the value provided 
by the model. Considering all the 16 results of tables III and IV , obtained using three 
free parameters: k b ,k c and 1 / z 0 we find 5 rms = 30%. It is important to remark that we 
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are fitting two different properties, the mass and the decay constant, of four resonances of 
two different hadrons using only three parameters in the model. So, the error of 30 % is 
reasonable, considering the simplicity of the model and the complexity of the strong coupling 
physics involved in the heavy meson structure. 


IV. FINAL COMMENTS 

The experimental data for the first four S-wave resonances of charmonium and bottomo- 
nium lead to decay constants that exhibit a clear decrease with the radial excitation level. 
In this article we proposed a consistent AdS/QCD approach where the decay constants are 
calculated from two point operator products at a finite value of the radial coordinate of 
anti-de Sitter space. The model was formulated with three independent parameters and was 
applied to the calculation of sixteen independent observables: eight masses and eight decay 
constants. 

A consistent AdS/QCD method for finding decay constants that decrease with radial 
excitation was lacking in the literature. This property is reproduced in our model, in contrast 
to the original soft wall, providing a more realistic description of heavy vector mesons. 
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los Andes. 
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